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We study the optimized version of the multiple invasion percolation model. Some topo- 
logical aspects as the behavior of the acceptance profile, coordination number and vertex 
type abundance were investigated and compared to those of the ordinary invasion. Our 
results indicate that the clusters show a very high degree of connectivity, spoiling the 
usual nodes-links-blobs geometrical picture. 
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1. Introduction 

When a nonviscous liquid is injected into a porous medium already filled with a 
viscous fluid the system can be found in two different regimes: one where the 
dominant forces are of capillary nature and another where the viscous forces are 
predominant. The theoretical descrintion of such system is based on two models: 
diff'usion-limited aggregation (DLA) a and invasion percolation □. DLA describes 
the fast displacement situation when the viscous forces are dominant. The invasion 
percolation model is applied to a slow fluid flow, where the leading forces are of 
capillary nature. The fluid displacement follows minimum resistance paths: the 
smaller pores are filled or invaded first. 

The original invasion percolation model have been modified in order to bring it 
as close as possible to real world. So, the action of an external gravitational field 
and the flux with a privileged direction LrLl were incorporated into the model. 
Recently □ the multiple invasion percolation model was proposed. In this model a 
certain number of lattice sites can be simultaneously invaded. It is always important 
to have in mind that invasion percolation is a dynamical model, that is, it intends 
to explain not only static properties like the clusters fractal dimensions but also 
how these clusters evolve in time. The lack of a meaning on associating or trying to 
make a correspondence between one step growth interval (of the algorithm) and one 
time interval separating two sucessives pore's invasions (of the real world) is what 
we were concerned by the time we formulated the multiple invasion percolation 
model a. In our models, one real time interval corresponds exactly to one growth 
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step and it is, in this sense, more realistic. When water penetrates layered soil, in 
some cases, a fingering formation is detected. The multiple invasion percolation not 
only describes successfully the fingers dynamics but also led us to a more ambitious 
program - to obtain the experimental pores size distribution of the soil Briefly, 
the multiple invasion percolation model is useful not only to explain experimental 
results but also to predict them. 

There are two kinds of multiple invasion percolation: the perimeter and opti- 
mized models. In the first model the cluster growth is controlled by the flux through 
the perimeter. The optimized model is governed by a scaling relation between the 
mass and the gyration radius of the cluster. Both types of invasion were studied 
in their site versions. For the perimeter model, the acceptance profile, mean coor- 
dination number and abundance of vertice type were investigated. In this paper 
we discuss the acceptance profile and some topological properties of the optimized 
model clusters. 

2. The Optimized Model 

We briefly recall the growth mechanism of the optimized modefi. It was devised 
to obey exactly the scaling 

M^iR,)^ (2.1) 

or as near it as possible {Rg is the gyration radius and D is a real positive external 
parameter that can be tuned). 

Basically we use the following strategy: at each growing step we build a list 
containing all the cluster perimeter sites that can be invaded and we ask for the 
number of sites that should be invaded in order that equation (2.1) is verified as 
osely as possible. When D G [1.89,2] these proceedings build a fractal object 
which is extremely stabilized (in the sense that in any stage or size the scaling 
is perfectly obeyed and not only in the asymptotic limit). In this interval the 
parameter D coincides with the fractal dimension of the cluster. Reference □ stress 
that although the ordinary invasion percolation and the optimized model (D = 1.89) 
have the same fractal dimension the clusters are not the same. While in the ordinary 
invasion percolation only one site is invaded at each step, in the optimized model 
several sites can be invaded at same time. The result is that the optimized model 
generates more massive clusters than those of the ordinary invasion percolation. 

In the region [0,1.89] the system is frustrated from below, that is, the scaling 
relation recommends an invasion of less than one site which is forbidden by the 
algorithm (at least one site must be invaded). In the region D > 2, the system 
is frustrated from above and a very beautiful burst phenomenon takes place u. In 
both cases, D does not coincide with the real fractal dimension. 

3. Numerical Simulations 



The Optimized Model of Multiple Invasion Percolation 3 



0.7- 



0.6- 



AA O OA OOO O A OAOOTAO O 

OfcAVOWWA* ■Aa»*AA AAAAAAA^»AAAO*AAAAO'»AO&A' BOA& 



S 0.5 H 

3^ 

CD 

0.4 H 



CD 

o 

-I—' 

Q. 
CD 
O 
O 
< 



0.3- 



0.2- 



0.1 - 



0.0- 



D = 1.8 



0.8- 
0.7- 
0.6- 
0.5- 
0.4- 
0.3- 
0.2- 
0.1 ■ 
0.0- 



L = 201 



D = 2.00 



D = 1.95 



— L = 51 
— A— L= 101 
— °— L= 151 
— L = 201 



\ 



• I A ▼ 

\ \ \ \ 



0.6 



1.0 



0.0 



I 

0.2 



I ' I 

0.4 0.6 

r 



0.8 



1.0 



Figure 1: Some acceptance profiles of the optimized model for fixed D and various 
lattice sizes. In the inset, we show the behaviour of a(r) for L = 201 and different 
values of D. 



In this paper we have studied the optimized model only in the region [1.89,2]. 
To analyze growth mechanisms it is useful to define i the acceptance profile a(r) 
which is the ratio between the number of random numbers in the interval [ r, 
r+dr ] accepted into the cluster and the total number of random numbers in that 
range. In the limit of an infinite lattice, the acceptance profile of the ordinary 
invasion percolation tends to a step function with the discontinuity located at the 
critical ordinary percolation threshold Pc = 0.5928 O. For the optimized model we 
determined the acceptance profile as a funtion of D. Fig. 1 shows the case D = 1.89 
for many lattice sizes. It does not seems to approach a step function, exhibiting 
a persistent tail near the inflexion point. This happens because, in order to obey 
the scaling rule, sites with larger random numbers are invaded. As a consequence, 
we find that the inflexion point position is bigger than its value in the ordinary 
invasion. The inset of the Fig. 1 shows that the inflexion point position increases 
with D. 

The sites of the multiple invasion percolation cluster can be classified according 
to their number of first nearest neighbors occupied sites. Let A''^ be the number 
of sites surrounded by k ( k — 1,..,4 for the square lattice) occupied sites. We 
determined the vertice concentrations of kind fc, i.e., ni- = -r-^^ — , for sizes L ~ 

51,101,151 and 201. The result was extrapolated using the BST algorithm 113 
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Table 1: The vertice concentrations Uk = ^ 'I, and the coordination number Z 

of the optimized model for many values of D. The data were extrapolated using 
the BST algorithm. 



and it is presented in the Table 1. The BST is a useful algorithm to extrapolate 
physical quantities that converge obeying a power law F{L) = F{oo) + AL^^ . It 
allows a reliable determination of critical parameters in the thermodynamic limit 
and its versatility becomes more pronounced if there are only very short sequences 
available. Our estimated coordination number for the optimized model {D = 1.89) 
is Z = 2.71 in contrast with Z = 2.51 for the ordinary invasion percolation i. In the 
table, we see that, as we increase D, the cluster becomes more compact favouring 
vertices of kind 4 in detriment of those of kind 1 which are practically extinct. 

The topological properties of the percolating cluster have been investigated for 
more than two decades. It is well known that the backbone is an important structure 
of the cluster. It is formed by the union of all self-avoiding walks connecting two 
points Pi and P2 of the lattice. This means that if we pass a current between 
these points, the backbone is identified as the set of sites that carries current. The 
backbone study has possible applications on the condutivity of random systems Il3 
and on fluids flowing in porous media 113. The backbone of tha^critical percolation 
cluster is non-Euclidean with fractal dimension 1.647 ± 0.004 llJ. 

Many pictures of the percolating cluster backbone were proposed 000. The 
most promising model was introduced by Stanley £3 and is know as the nodes-links- 
blobs model. In this model the backbone consists of a network of nodes connected by 
one-dimensional links which are often separated by multiconnected pieces or blobs 
of all length scales lia. Thus, the backbone may be viewed as a topologically linear 
string of blobs of all possible sizes. In general, there are sites that, when removed, 
split the spanning cluster into pieces. These sites are termed red sites. In the nodes- 
links-blobs model the red sites are also called blobs of size one E3. At criticality, the 
number of red sites scales as a power law with the exponent calculated exactly by 
ConigHo e1 i = 0.75. 

We have already studied the backbone and elastic backbone structures of the 
multiple invasion percolation model lEI. To determine these structure we employed 
the burning algorithm E3. We chose this algorithm because, beyond the backbone 
and the elastic backbone, it also permits the determination of the red sites number. 
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Figure 2: A typical cluster of the optimized model for D = 1.89 and L = 201. The 
cluster has 14650 sites. 

the minimum path and the loops number. 

For the optimized model, the scaling (equation (2.1)) is also perfectly obeyed by 
the backbone. At D = 1.89 we got the fractal dimension Dp = 1.74 ± 0.01. This 
means that although the optimized model at D = 1.89 and the ordinary invasion 
have the same fractal dimension, they arc intrinsicaly different since their backbones 
are not the same. The number of red sites A^^ is very small and random. It does not 
seem to obey any power law. The red sites number is so small that the probability 
of disconnecting the cluster by removing any site randomly is pratically zero. The 
optimized algorithm destroys the red sites, increasing the cluster connectivity. The 
disappearence of the red sites indicates that the nodes-links-blobs model may not 
be useful to describe the backbones of the optimized model clusters. Indeed, the 
cluster of the optimized model has a very different form from that of the ordinary 
invasion percolation as can be seen in the Fig. 2. 

4. Conclusions 

We studied the optimized model of multiple invasion percolation by comparing 

its topological properties with those of the ordinary invasion percolation. The high 
connectivity of the clusters produced by the optimized model together with the 
almost disappearance of the red sites, led us to conclude that the nodes-links-blobs 
model is not well suited to describe the backbones of the optimized model. 

5. Acknowledgements 



6 R. A. Zara & R. N. Onody 



Wc acknowledge CNPq (Conselho Nacional de Desenvolvimento Cientifico e Tec- 
nologico) and FAPESP ( Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
) for the financial support. 

References 

1. T. A. Wittcn and L. M. Sanders, Phys. Rev. B 27,5686 (1983). 

2. D. Wilkinson and J. F. Willomscn, J. Phys. A 16, 3365 (1983). 

3. D. Wilkinson, Phys. Rev. B 30, 520 (1984). 

4. A. Birovljev, L. Furuberg, T. Feder, T. Jossang, K. J. Mal0ay and A. Aharony, Phys. 
Rev. Lett. 67 584 (1991). 

5. P. Moakin, J. Foder, V. Frette and T. Jossang, Phys. Rev. A 46, 3357 (1992). 

6. R. N. Onody, A. N. D. Posadas and S. Crostana, J. Appl. Phys. 78, 2970 (1995). 

7. R. N. Onody, Int. ,J. Mod. Phys. C 6, 77 (1995). 

8. R. N. Onody and R. A. Zara , Physica A 231, 375 (1996). 

9. B. B. Mandelbrot, T/ie Fractal Geometry of Nature Freeman, New York (1982). 

10. R. Bulirsch and J. Stoer, Numor. Math. 6 (1964) 413. 

11. D. StaufFer and A. Aharony in Introduction to Percolation Theory Taylor and Fran- 
cis, London (1992). 

12. P. G. de Gennes, J. Physique Lett. 37, LI (1976). 

13. H. E. Stanley and A. Coniglio, Phys. Rev. B 29, 522 (1984) 

14. P. Grassberger, J. Phys. A 25, 5475 (1992). 

15. A. Skal and B. L Shoklovskii, Fiz.Tekh. Poluprovodn 8 (1974) 1586. [Sov. Phys. Semi- 
cond. 8 1029 (1975) 

16. Y. Gefen, A. Aharony, B. B. Mandelbrot and S. Kirkpatrick, Phys. Rev. Lett. 47, 
1771 (1981). 

17. H. E. Stanley J. Phys. A 10, L211 (1977) 

18. H. J. Herrmann and H. E. Stanley, Phys. Rev. Lett. 53,1121 (1984). 

19. A. Coniglio Phys.Rev. Lett. 62, 3054 (1989) 

20. R. N. Onody and R. A. Zara , Phys. Rev. E 56, 2548 (1997). 

21. H. J. Herrmann, D. C. Hong and H. E. Stanley, J. Phys. A 17, L261 (1984). 



